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ABSTRACT 

We explore the relative strength of the narrow emission lines in an SDSS based sam- 
ple of broad Ha selected AGN, defined in paper I. We find a decrease in the narrow to 
broad Ha luminosity (LbHa) ratio with increasing LbHa, such that both L([0 III] A5007) 
and L(narrow Ha) scale as °< L^ a for 10 40 < L b H« < 10 45 erg s _1 . Following our ear- 
lier result that LbHa x Lboi, this trend indicates that the relative narrow line luminosity de- 
creases with increasing Lboi- We derive L bo i/10 43 erg s -1 = 4000(L([O Hl])/10 43 erg s -1 ) 139 . 
This implies that the bolometric correction factor, L DO i/L([0 III]), decreases from 3000 at 
Lboi = 10 461 erg s _1 to 300 at Lboi = 10 42 5 erg s _1 . At low luminosity, the narrow compo- 
nent dominates the observed Ha profile, and most type 1 AGN appear as intermediate type 
AGN. Partial obscuration or extinction cannot explain the dominance of intermediate type 
AGN at low luminosity, and the most likely mechanism is a decrease in the narrow line region 
covering factor with increasing Lboi- Deviations from the above trend occur in objects with 
L/Ledd ^ 10~ 26 , probably due to the transition to LINERs with suppressed [O III] emission, 
and in objects with Mbh > 10 8 5 M Q , probably due to the dominance of radio loud AGN, and 
associated enhanced [O III] emission. 
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1 INTRODUCTION 

The circumnuclear gas located on 1 - 1 000 pc scale from active 
galactic nuclei (AGN) is an important constituent in several intrigu- 
ing, but poorly known processes. It is the source of AGN fuel, and 
it absorbs AGN energy and momentum via winds, radiation and/or 
jets. Understanding the properties of this gas is crucial to constrain 
theories of these processes, and their possible connection to star 
formation in the host bulge, which is implied from the central black 
hole - bulge relations (Magorrian et al. 1998, Ferrarese & Merritt 
2000). 

Some of the properties of this gas have been inferred by ex- 
amining how photons from the central source are reprocessed into 
secondary radiation. The AGN radiation is mainly converted either 
into thermal IR emission from dust grains embedded in the gas, or 
to permitted and forbidden emission lines. These emission lines are 
observed with widths of ~ 300 km s , typical of the galaxy poten- 
tial. The region that emits the lines is termed the narrow line region 
(NLR), as opposed to the broad line region (BLR, ~ 3 000 km s" 1 ), 
which is located on smaller scales, within the black hole gravita- 
tional sphere of influence. The NLR is used to constrain the radial 
distribution of the circumnuclear gas (Ferguson et al. 1997), and 
to infer the AGN bolometric luminosity Lboi when other measures 
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are not available (e.g. Kauffmann & Heckman 2009; Schawinski et 
al. 2010). This circumnuclear gas may also block our line of sight 
to the inner ionizing source and BLR, thus playing a role in the 
unobscured (type 1) / obscured (type 2) AGN dichotomy. 

How does the distribution of the circumnuclear gas depend 
on Lboi an d black hole mass Mbh? There is some evidence that 
the covering factor CF of the IR emitting gas decreases with Lb i 
(Maiolino et al. 2007, Treister et al. 2008), as does the type 2 / type 
1 ratio (Hasinger 2008). The dependence of the NLR covering fac- 
tor, CFnlr, on luminosity has not been directly explored, though 
there may be some indications that CFnlr is higher at low AGN 
luminosity (Ludwig et al. 2012). Also, at high Lbol> where the host 
contribution to the continuum is negligible, a decrease in [O III] 
equivalent width with luminosity has been observed (Sulentic et 
al. 2004, Netzer et al. 2006, Ludwig et al. 2009), which may also 
imply a luminosity-dependent CFnlr- 

Moreover, in the classical type 1 / type 2 AGN division, the 
broad lines either dominate the permitted lines or are not observed 
at all. This division is unsatisfactory, as the broad and narrow com- 
ponents of the permitted lines show a large range in flux ratio be- 
tween different objects (e.g. fig. 1 of Stern & Laor 2012, hereafter 
Paper I). Following Osterbrock (1977), 'intermediate types' were 
added to the scheme - type 1.2s, 1.5s and 1.8s, with larger num- 
bers corresponding to higher narrow to broad flux ratios. If in type 
1 AGN we have a direct view of the BLR, and in type 2 AGN the 
BLR is obscured, what then are the intermediate types? A simple 
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interpolation of the above scheme may imply that the BLR is either 
partially obscured or undergoes extinction by optically thin dust. 
Partial obscuration is clearly a factor in AGN, as can be seen in 
spectropolarimetric surveys (e.g. Tran 2003) that show objects in 
which a 'hidden' BLR is revealed in the polarized spectrum. Evi- 
dence for dust reddening of the BLR has also been found in inter- 
mediate type samples (e.g. Goodrich et al. 1995). However, Trippe 
et al. (2010) find that the majority of intermediate types in their 
sample are inconsistent with reddening of the BLR, which suggests 
another mechanism may be important in driving the intermediate 
type phenomenon. An alternative mechanism is that intermediate 
type AGN are not obscured, like 'standard' type Is, but have a large 
CFnlr/CFblr ratio. As we show below, this mechanism appears 
to dominate in low luminosity AGN. 

In this paper, we use the type 1 sample described in Paper I. 
The sample is based on the detection of broad Ha emission, which 
in some objects is relatively weak compared to the narrow Ha 
emission, classifying them as intermediate type AGN. The sam- 
ple spans 10 6 < M BH < 10 9 5 M Q and 10 42 < L bol < 10 46 erg s _1 . 
We derive the NLR luminosity and fraction of intermediate types, 
for different L bo i and M BH . Using these relations between the NLR 
and AGN properties, we explore various mechanisms proposed to 
explain the intermediate type phenomenon. We also calibrate the 
luminosity-dependent ibol/^[Om] ratio, for use of future studies 
based on type 2 AGN. 

The paper is organized as follows. In §2, we give a summary 
of the sample selection and emission line measurements described 
in Paper I, emphasizing the measurements of [O III] and narrow 
Ha, used here. In §3 we present our main result, the increase in 
relative NLR luminosity and in the fraction of intermediate type 
AGN with decreasing luminosity. The relative NLR luminosity is 
also compared to L/L^dd an d A^bh- In §4 we compare our results 
to earlier studies. The possible physical mechanisms behind the ob- 
served trend are examined in §5. We discuss the implications in §6, 
and provide conclusions in §7. 

Throughout the paper, we assume a FRW cosmology with fl m 
= 0.3, A = 0.7 and Hq = 70 km s _1 Mpc" 1 . 



2 THE SAMPLE 

The construction of the broad Ha selected sample used in this pa- 
per is described in Paper I. In §2.1 we provide a brief summary, 
emphasizing the details most relevant to this work. In §2.2, we dis- 
cuss a change in the host subtraction procedure. In §2.3 and §2.4 
we elaborate on the measurement of the [O III] luminosity Ljo m] 
and the narrow Ha luminosity L n H«, which are the basis of the 
analysis in this paper, and were not used previously. 

2.1 The Tl sample creation 

The type 1 (Tl) sample is selected from the 7 th data release of 
the Sloan Digital Sky Survey (SDSS DR7; Abazajian et al. 2009). 
The SDSS obtained imaging of a quarter of the sky in five bands 
(ugriz; Fukugita et al. 1996) to a 95% r band completeness limit 
of 22.2 mag. Objects are selected for spectroscopy mainly due to 
their non-stellar colors (Richards et al. 2002), or extended mor- 
phology (Strauss et al. 2002). The spectrographs cover the wave- 
length range 3800A-9200A at a resolution of ~ 150 km s" 1 , and 
are flux-calibrated by matching the spectra of simultaneously ob- 
served standard stars to their PSF magnitude (Adelman-McCarfhy 
et al. 2008). 



We use SDSS spectra which are classified as non-stellar and 
have 0.005 < z < 0.31. To ensure a reliable decomposition of the 
broad and narrow components of Ha, we use only spectra with 
S/N > 10 and sufficient good spectral pixels in the vicinity of 
Ha, as detailed in Paper I. These requirements are fulfilled by 
232 837 of the 1.6 million spectra in DR7. The spectra are cor- 
rected for foreground dust, using the maps of Schlegel et al. (1998) 
and the extinction law of Cardelli et al. (1989). We then subtract 
the host, as detailed below. Then, we subtract a featureless contin- 
uum, derived by interpolating the mean continuum level at 6 125 A- 
6250A and 6880A-7000A. The residual flux at 6250A-6880A 
(±14,000 km s _1 from Ha) is then summed, excluding regions 
±690 km s- 1 from the [O I] 116300,6363, [N II] 116548,6583, 
[S II] 116716,6730 and Ha narrow emission lines. We find 6986 
objects with significant residual flux, which is potential broad Ha 
emission. 

On the objects with significant residual near Ha, we fit the 
profiles of the broad and narrow Ha, [O III] 15007, and the [O I], 
[N II] and [S II] doublets mentioned above. Narrow lines are fit 
using 4 th -order Gauss-Hermite functions (GHs; van der Marel & 
Franx 1993) and an up to 10 th -order GH is used for the broad Ha 
profile. The following criteria are applied to the broad Ha fit, in or- 
der to exclude objects in which the residual flux is not clearly BLR 
emission: the FWHM (Av) of the fit is required to be in the range 
1000 - 25000 kms" 1 ; the total flux of the fit, and its flux den- 
sity at the line centre, are required to be significant. The Tl sample 
consists of the 3 579 objects which passed these criteria. The broad 
Ha luminosity (L b ea) an d Av of the 3 410 Tl objects used here 
(see below) are listed in Table 1. The selection effects implied by 
our selection criteria are detailed in Paper I. 

We supplement the optical SDSS spectra in the Tl sample by 
photometric measurements in the UV and X-ray, as described be- 
low. 



2.1.1 L uv 

The GALEX mission (Martin et al. 2005, Morrissey et al. 2007) 
observed 2/3 of the sky in the FUV (effective wavelength 1528A) 
and NUV (2271 A) bands. We search the GALEX GR6 catalog for 
objects within 5" of the Tl objects. Of the 89% of the Tl objects 
observed by GALEX, 93% have detections in both bands. We use 
UV fluxes derived from the 6"-radius 'aperture 4', corrected for 
the PSF and for foreground dust assuming Anuv/^b-v = 8-2 and 
Afuv/£b-v = 8.24 (Wyder et al. 2007). In case of multiple detec- 
tions per object, we co-add the observations, weighted by exposure 
time. For an object that was not detected, we assign an upper limit 
equal to the mean upper limit of the relevant GALEX survey (Mor- 
rissey et al. 2007), scaled to the exposure time at the field of the ob- 
ject. We derive Luv = vL v (1528A), the luminosity of the shortest 
restframe wavelength available in all detected objects, by a power 
law interpolation between the two UV bands. The Luv values are 
listed in Table 1. 



2.1.2 L x 

The ROSAT All-Sky Survey (Voges et al. 1999, 2000) covers the 
entire celestial sphere in the 0.1— 2.4 keV range, with positional 
uncertainties of 10" - 30". We find matches to 1 561 (43%) of the 
Tl sample within 50" of each Tl object, expecting 7 objects to 
have false matches. Only 0.4% of the Tl objects were targeted for 
spectroscopy by the SDSS solely for being near a ROSAT source. 
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Therefore, we do not expect our sample to be biased towards X- 
ray bright AGN. The ROSAT count rates are converted to L x = 
vL v (2 keV) assuming a power-law X-ray spectrum with cc x = 1.5 
(Laor et al. 1994, Schartel et al. 1996), the z of the optical match, 
and the Galactic Nn measurements of Stark et al. (1992). Upper 
limits in Tl objects without ROSAT detections are set to the typical 
limiting sensitivity of 1CT 12 ' 5 erg s _1 cm~ 2 . The Lx values are 
listed in Table 1. 



2.2 Host subtraction 

In Paper I, we fit three galaxy eigenspectra (ESa) from Yip et al. 
2004, and aLj «1 -15 power law representing the AGN contin- 
uum, to wavelength regions devoid of strong emission lines (Paper 
I, appendix A2). We then account for the non-Balmer absorption 
lines near Ha by subtracting the three fit ESa. However, since the 
Yip et al. ESa include emission lines which we interpolate over 
prior to the subtraction, the B aimer absorption features are not ac- 
counted for. This prescription sufficed for modeling the broad Ha 
component, as the observed equivalent width (EW) in the Tl sam- 
ple is typically ~ 100A, and > 15A in 99% of the objects. In con- 
trast, the Ha stellar absorption EW is typically lk for old stars (see 
fig. 1 in Hao et al. 2005), and < 5A for a young stellar population 
(Groves et al. 2012). 

However, the apparent flux of the weaker narrow Ha might 
be more strongly affected by stellar Ha absorption. Therefore here, 
instead of a simple interpolation, we replace the emission lines in 
ESI with absorption features of old stars from the Hao et al. (2005) 
ES 1 . The latter was derived from a principle component analysis of 
galaxies without emission lines. The fact that the continua in both 
the Yip et al. ESI and the Hao et al. ESI are dominated by old 
stars justifies the replacement. Technical details of this procedure 
are given in appendix A. 

Our fit does not account for the possible Balmer absorption 
from young stars. Their possible effect on L n H« is addressed below. 

We note that the Ha absorption feature in Hao ESI has fj ~ 
400 km s _1 , compared to the typical a ~ 150 km s _1 of other ab- 
sorption features, indicating the intrinsic width dominates the width 
of the feature. Therefore, broadening the ES during the fit will have 
a negligible effect on the Ha absorption feature in all but the most 
massive galaxies. Hence, we avoid this additional complexity in the 
fitting algorithm. 

How does this adjustment of the Yip et al. ESI affect the mea- 
sured broad Ha? In 98% of the objects, the change in L^ea an d 
Av, compared to the Paper I result, is < 0. 1 dex. This negligible 
change is expected due to the large broad Ha EW. From inspec- 
tion, 58 of the remaining 2% (74 objects) have reasonable fits. The 
remaining 16 objects were not further processed due to their rela- 
tively small number, and they were removed from the sample. Four 
objects were fit with Av < 1000 km s~ 1 , below our selection thresh- 
old, and therefore removed as well. 



2.3 The [O in] measurements 

Since the [O III] line is used extensively in this paper, we use only 
objects which have sufficient good pixels in its vicinity to ensure a 
reliable fit, as detailed in appendix A4.1 of Paper I. The 149 objects 
that fail this criterion are removed from the sample. The final Tl 
sample size is 3 410 objects. 

The profile of the [O III] A5007 emission line can be blended 
with the adjacent [Fe II] complex and with the red wing of the broad 



SDSS Name 


L btia 


Av 


^[O III] 


L nHcc 


L uv 


w 


Notes 


J000202.95- 103037.9 


41.9 


2310 


41.4 


41.5 


43.8 


42.6 




J000410.80-104527.2 


42.6 


1360 


41.5 


41.9 


44.6 


43.2 




J00061 1.55+145357.2 


42.1 


3320 


40.6 


41.0 


44.0 


42.6 




J000614.36-010847.2 


41.6 


3910 


40.9 


41.1 


43.2 


41.9 


-,-,-,u 


J00065776+152550.0 


41.5 


3020 


40.8 


40.7 


43.0 


41.7 


-,-,-,u 



Table 1. The broad and narrow line measurements and photometric vL v lu- 
minosities of the Tl sample. All luminosities are in log erg s _1 , while Av 
is in km s _1 . The Luv and L\ are measured at rest wavelengths of 1528A 
and 2 keV. Notes for L[o m] , LiHa > Luv an d Lx are separated by commas, 
and coded as follows: '-': no note; 'U': upper limit; 'E': narrow line mea- 
surement relative error is > 30%; 'O': narrow Ha is [O III]-like (§2.4.1); 
'A': Lniia measurement is affected by stellar absorption (§2.4.2); 'N': not 
observed by GALEX. The electronic version includes all 3 410 objects. 

H/3 profile. Therefore, we fit the wavelength region 4967A-5250A 
with three components: an iron template derived from observations 
of I Zw 1 (kindly provided by T. Boroson), a GH for the [O III] line, 
and a broken power law for the underlying continuum and H/3 red 
wing. The measured Ljo m] of the Tl objects are listed in Table 1. 
In 99% of the objects, the relative error on the [O III] flux F[q m] is 
< 30%. In 0.5% of the objects [O III] is not detected, and the quoted 
values are upper limits. The derivation of the error and upper limits 
is given in appendix Bl. 

We compare our results with the Ljom] measurements of 
SDSS quasars in Shen et al. (2011), which model [O III] and H/3 
using multiple Gaussians. The blending of the [O III] profile with 
the [Fe II] multiplets and broad H/3 is pronounced mainly in lu- 
minous AGN (see fig. 18 in Paper I). Therefore, the discrepancy 
between different deblending algorithms in Tl objects that appear 
in the SDSS quasar catalog should be an upper limit on the discrep- 
ancy in the whole Tl sample. In the 419 overlapping objects, the 
mean ratio of our L[o m] to the Shen et al. F[o m] is 0.04 dex, with 
a dispersion of 0.06 dex. 



2.4 The narrow Ha measurements 

The narrow and broad components of Ha are fit simultaneously 
with the adjacent [O I], [N II] and [S II] forbidden lines, keeping the 
width and higher GH coefficients (A3, of all narrow lines equal, 
and using a 4 th -order GH for the broad Ha. We perform three fit 
attempts, in which and of the narrow lines are either zero (a 
pure Gaussian), same as in the [O III] profile, or can vary between 
these two values. We use the result with the lowest y}, and refine 
the fit by adding GH coefficients to the broad component. Further 
details of the algorithm are given in appendix A4.2 of Paper I. In 
appendix B2 here, we show that the three different assumptions on 
hj, and A4, and the different number of GH coefficients used for the 
broad component, have a negligible effect on the measured narrow 
Ha flux, F Ma . 

The measured L n H« of the Tl objects are listed in Table 1. Two 
sources of uncertainty in L n ua are described in the two following 
subsections. Additionally, objects in which only an upper limit on 
F n n a can be measured (2% of the sample), and objects with > 30% 
relative error (3%), are described in appendix B2. 



2.4. 1 [O III ]-like narrow Ha 

In 15% of the sample, the fit yielded FWHM(nHa) ^ 1.5x 
FWHM([0 III]). These objects have non- or barely-detectable nar- 
row lines near Ha, and there is no clear transition between the 
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Figure 1. Luminosity vs. FWHM distribution of the total (broad + narrow) Ha emission line in the 3 410 objects of the Tl sample. Yellow triangles indicate 
Tl objects in which Av t01a i is overestimated due to stellar Ha absorption. With decreasing Ln a , Av to tai drops sharply, reaching typical NLR widths of 
~ 300 km s~ 1 . Almost all log L\\ a > 43 objects have Av tola i > 1 000 km s~ ' , while most objects with log Lira < 41 .5 have Av tota i < 1 000 km s~ 1 . Furthermore, 
most low luminosity high Av tota i objects are upper limits. The dominance of low Av tota i profiles at low luminosity is due to an increase in the L„Ha /^bHa ratio 
with decreasing luminosity. Therefore, selecting broad line AGN using a Av lota i > 1 000 km s criterion (e.g. Vanden Berk et al. 2006, Schneider et al. 2010) 
will be highly incomplete at low AGN luminosities. 



broad and narrow components of Ha. In such cases, we fit the nar- 
row Ha using the FWHM, hi and h\ of [O III]. As F n e« is less 
certain in these objects, they are marked by different symbols in 
figures where F n ua ls used. 

It is possible that with higher spectral resolution data the de- 
marcation between the BLR and NLR in these objects would be 
more clear, and one could test if the proposed fits remain consis- 
tent. 



2.4.2 Stellar absorption 

The Ha absorption feature of young stars, which has EW< 5A 
(Groves et al. 2012), is not modeled by our fitting algorithm. There- 
fore, when the ratio of AGN to host luminosity is low, the mea- 
sured F„Ha and Ha flux density can be underestimated due to the 
underlying stellar absorption. We note that the typical absorption 
FWHM is 1 300 km s -1 (Yip et al. 2004), compared to the typical 
FWHM(nHa) of 300 km s~ 1 . Therefore, the error in is likely 
< 3 A x L A (host). 

We derive the galaxy continuum luminosity density in the 
Tl objects by subtracting the AGN contribution, using EW(bHa) 
= 570A (Paper I). In 156 (5%) of the Tl objects, L nHa < 3A x 
L^(host), implying a potentially strong error on L n Ha due to stellar 
absorption. The L n ^ a in these objects are marked as lower limits. 
Also, in these objects we mark the FWHM of the total Ha profile 
(Avtotai) as an upper limit, since the measured flux density of the 
total Ha is a lower limit. 



found that L^ a i = 130* 24 A)Ha- Thus, we use LbHa as a proxy for 
Lbol> and LbHa and Av to derive Mbh (see eq. 2 there). We begin 
by analyzing the total Ha line, i.e. the sum of the broad and nar- 
row components, as it is less sensitive to errors in the deblending 
procedure. We then proceed to examining L nHa and m] directly. 

3.1 The total Ha 

Figure 1 shows the distribution of the Tl objects in the total 
Ha FWHM (Av tota i) vs. total Ha luminosity (Lh«) plane. At 
Lhc > 3 x 10 43 erg s _1 , practically all Tl objects have Av tota i > 
1000 km s , indicating the broad component dominates the 
Ha profile. With decreasing luminosity Av tota i decreases, such 
that 85% of the Tls with L Ha < 10 42 erg s -1 have Av tota i < 
1 000 km s" 1 . At L Ha < 10 41 erg s" 1 , 73% of the Tl objects have 
Av tota i < 500 km s . Widths below 500 km s are typical of 
the NLR, indicating that at low luminosities the NLR dominates 
the Ha profile, or equivalently, most AGN are intermediate types. 
Some of the objects with low L Ha do have high Av tota ], but almost 
all of them have a stellar absorption feature with strength compa- 
rable to the narrow emission line (§2.4.2), indicating their Av tota i 
are upper limits, and reinforcing the observed trend of Av tota i with 

In previous works, broad line AGN were selected based on 
optical spectra using a Av tota i > 1000 km s criterion (Vanden 
Berk et al. 2006, Schneider et al. 2010). Figure 1 implies that this 
criterion is adequate for quasars, but is highly incomplete at low 
AGN luminosities. 



3 THE RELATIVE NLR LUMINOSITY 



Our goal is to understand the dependence of the fraction of inter- 
mediate type AGN on AGN characteristics. In Paper I, we have 
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Figure 2. Mean 6280A-6800A spectra in bins of LbHa and Av, overplotted on the LbHa vs. Av distribution of the Tl sample. In each bounding box (dotted lines) we plot the mean observed spectrum, before host 
subtraction, of all objects with LbHa and Av in the range delimited by the box. The y-axis of each spectrum is in erg s A , and extends from the continuum level to the max flux density of Ha, so Av lota i (Fig. 1) is 
measured at the middle of each bounding box. The prominent features seen in the spectra are the broad and narrow Ha, and the [O I] AA6300, 6363, [N II] AA6548, 6583 and [S II] AA6716, 6731 low ionization 



forbidden lines. Note the increasing prominence of the narrow Ha with decreasing luminosity, as indicated by Fig. 1. At LbHa ^ 10 



the mean Ha profiles are typical of Seyfert 1.0s, at LbHa ~ 10 er g s 



the mean Ha profiles are typical of Seyfert 1.5s, while at LbHa ~ 10 40 — 10 41 erg s _1 type 1 AGN tend to be Seyfert 1.8s. The binning by LbHa and Av is equivalent to binning by Mbh and L/Lem- The Mbh increases 
upward and rightward, while L/L E(M increases upward and leftward. At high L b Ha and low Av, or at high L/Le<m, there is no clear transition between the narrow and broad components of the Ha profiles. In such 
objects, we use the profile of [O III] A50O7 as a template for the narrow Ha (§2.4.1). With decreasing L/Lem, the low ionization forbidden lines become stronger relative to the narrow Ha, probably indicating a 
transition to LINER narrow line ratios (see §3.4.3). 
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In Figure 2, we plot the mean spectra of the Tl objects in bins 
of LbHa and Av, m the range 6280A-6800A. Each bin spans a fac- 
tor of ten in LbHa> and a factor of two in Av. The peripheral bins 
are not evenly populated. We therefore show in the background the 
distribution of the objects in the LbHa vs. Av plane. Each mean 
spectrum is derived from all objects inside its bounding box. Mean 
spectra are calculated by geometrically averaging luminosity den- 
sities of spectrum pixels with the same restframe wavelength X, 
rounded to 1CT 4 in logA. The spectra that enter this averaging are 
observed spectra, before host subtraction. The y-axis of each spec- 
trum extends from the continuum level to the maximum flux den- 
sity of Ha, so Av tota i (Fig. 1) corresponds to the width of Ha at the 
middle height of each bounding box. 

With decreasing luminosity, the narrow Ha component dom- 
inates the total Ha profile, as indicated by the Av tota ] analysis in 
Fig. 1. At LbHce > 10 43 erg s _1 (top row), the mean Ha profiles 
are dominated by the broad component, typical of Seyfert Is. At 
10 43 > LbHa ^ 10 42 er g s ~' (second row), the mean Ha pro- 
files indicate that most objects are Seyfert 1.5s, while at LbHa ~ 
10 40 — 10 41 erg s _1 (lowest row), the objects tend to be Seyfert 
1.8s. Intermediate type AGN dominate the population of broad line 
AGN at L bH a < 10 42 ergs" 1 . 

Note that a relatively stronger Ha stellar absorption feature 
at low luminosity, which is probably associated with the stronger 
observed stellar features, will only diminish the observed narrow 
Ha (see appendix A). Therefore, the true trend in narrow to broad 
ratio is likely stronger than observed in Fig. 2. 

In Fig. 2, Mbh increases upward and rightward, while L/Lgdd 
increases upward and leftward. The dependence of the NLR to BLR 
luminosity ratio on L/L Edd and M BH is examined below (§3.4). 
Note that at high L b Ha and low Av, or high L/L E( jd, the profiles 
show no clear transition between the narrow and broad components 
of Ha. In these objects, we fit the narrow Ha with a profile iden- 
tical to [O III] (§2.4.1). Also, with decreasing LbHa and increasing 
Av, or with decreasing L/L E ^, the low ionization forbidden lines 
become stronger relative to the narrow Ha. The relative ratios of 
narrow emission lines, and their dependence on L/L E( jd, M B h, L do i 
and Av, are the focus of a following paper. 

3.2 The narrow Ha 

We now utilize the L nHa measurements described in §2.4. In Fig- 
ure 3, the luminosity ratio of the narrow and broad Ha is plotted 
versus LbHa- Black dots mark objects with robust L n Ha measure- 
ments, while other colors and shapes indicate a possible bias in 
LnHa- A trend of increasing L n Ha/LbHa with decreasing luminos- 
ity is clearly seen, extending over a range of 10 4 in LbHa • The least- 
squares best fit power law of all Tl objects is 



LnHa; 42 = 0. 13 L^™ , O = 0.40 dex 



(1) 



where both luminosities are given in units of 10 42 erg s _1 , and a is 
the scatter of the measured log L n Ha around the relation. We treat 
L b Ha> which was used to select the Tl sample, as the independent 
variable. 

The [O IIl]-like narrow Ha are located at the high-LbHa end 
of the sample, with a L n Ha/LbHa ratio consistent with objects with 
robust L nHa measurements. Objects with strong stellar absorption 
are biased to lower L nH a than other Tl objects with the same LbHa> 
as expected if part of the narrow emission vanished in the Ha ab- 
sorption feature. Overall, the trend of L n Ha/LbHa with LbHa is in- 
dependent of the details of the Ha fit, as implied by Figs. 1 and 
2. 



Excluding objects with [O IIl]-like narrow Ha or strong stel- 
lar absorption, we get 



L„Ha; 42 = 0. 15 L^f™ 1 , a = 0.37 dex 



(2) 



similar to equation 1. The measurement error on L n Ha is < 0.1 dex 
(§2.4), therefore the noted dispersion is dominated by the intrinsic 
dispersion. 

Next, we evaluate the abundance of intermediate type AGN 
at different AGN luminosities. We define Seyfert 1.5s and 1.8s 
analytically, using the L n H a /LbHa ratio. The difference between 
Seyfert 1.8s and 1.9s is disregarded, as it is based on the detectabil- 
ity of the broad H/3. Objects with L nH a ~ 0.1 L b Ha> i- e - a roughly 
equal broad and narrow flux density for the typical width ratio of 
10 (= 3000/300), are considered as Seyfert 1.5s. Objects with 
LnHa ~ 0.3 LbHa (i- e - nux density ratio ~ 3) are considered as 
Seyfert 1.8s. 

Fig. 3 shows that at L bHa > 10 43 erg s -1 , only 20% of the 
objects are above the Seyfert 1.5 line, compared to 77% at L b Ha < 
10 42 ergs~'. At L b Ha < 10 41 erg s -1 , 52% of the objects are above 
the Seyfert 1.8 line. Intermediate types dominate the broad line 
AGN population at low luminosity, as already implied by Figs. 1 
and 2. 



3.3 The [O m] line 

We use the [O III] luminosity as another measure of the NLR lumi- 
nosity. Measuring the [O III] line is complementary to the narrow 
Ha measurement, as [O III] is less susceptible than Ha to con- 
tamination from H II regions powered by star formation in the host 
galaxy, but is more sensitive to ionization and density effects in the 
NLR. In Figure 4, we plot L[Om]/LbHa vs - LbHa in the Tl sam- 
ple. The least-squares best fit slope, using LbHa as the independent 
variable, is 



-[O ill] ; 42 



:0.16L 



0.72±0.01 
bHa;42 ' 



a = 0.36 dex 



(3) 



similar in slope to the L n Ha vs. L b Ha relation. The measurement 
error on L[o m ] is < 0.1 dex (§2.3), therefore the noted dispersion 
is dominated by the intrinsic dispersion. 

We add 20 234 z < 0.82 AGN from the SDSS quasar catalog 
(QCV, Schneider et al. 2010), in order to extend our luminosity 
range. We use the L bH| g measurements of Shen et al. (2011), con- 
verted to LbHa using LbHa/LbH/S = 3-2, the mean ratio for 419 
objects common to QCV and the Tl sample (the dispersion in 
LbHa/LbH/J is 0-1 dex). We also subtract 0.04 dex from the Shen 
et al. L[Om] (§2.3). To prevent confusion, Fig. 4 shows only the 
density contours of the QCV objects. These follow the same trend 
as the Tl objects, up to LbHa = 5 x 10 44 erg s _1 . The trend of in- 
creasing L[o in] /LbHa with decreasing L b Ha extends over 5 orders 
of magnitude in luminosity. 



As discussed in the introduction, Li 



[O in] 



is often used as 



a measure of AGN bolometric luminosity, when other measures 
are unavailable, particularly in obscured (type 2) AGN. Heckman 
et al. (2004) found Lb i/L[o ni] = 3 500, derived from the mean 
EW([0 III]) of z < 0.3 type 1 AGN. Here, we derive L bol /L[om] 
on the Tl sample, using the measured L[o m]/L b Ha> and the lumi- 
nosity independent Lb i/LbHa = 130 (Paper I). The derived relation 



Lbol 



= 4000^(1043 



Om] 



10 43 erg s _1 ^* ■ lir j erg s~ 

Derivation of the uncertainty is described below. 



(4) 
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Figure 3. The distribution of the NLR to BLR luminosity ratio in the Tl sample, vs. LbHa. Black dots mark Tl objects with robust L„Ha measurements, 
while other markers indicate the measured L„Ha could be biased. The first relevant uncertainty of the following list determines the shape and color of the 
marker: no clear BLR/NLR transition (red 'x', see §2.4.1); L n Ha is underestimated due to stellar absorption (yellow triangles, §2.4.2); L n Ha is an upper limit 
(blue arrows). The slope of the best fitting power law (dashed line) and the associated dispersion are noted. The tendency of increasing NLR/BLR ratio 
with decreasing luminosity implied by Figs. 1 and 2 is clearly seen, extending over a range of 10 4 in LbHa- Objects without a robust L„Ha measurement 
are consistent with the trend implied by objects with an accurate measurement. Objects near or above the Seyfert 1.5 line are intermediate type AGN. The 
abundance of intermediate types at low luminosities disfavors partial obscuration or variability as their origin. 
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Figure 4. The distribution of the [O III] A5007 to broad Ha luminosity ratio vs. LbHa, in the Tl and QCV samples. Tl objects are shown as black points, 
with upper limits on L[o mi marked by blue arrows. The slope of the best fitting power law (dashed line) and the associated dispersion are noted. The QCV 
distribution (measurements from Shen et al. 201 1) is shown as density contours, where at each 0.3 dec wide LbHa bin, 50% of the objects lie between the inner 
red dotted lines, while 80% lie between the outer yellow lines. The decrease in relative [O III] luminosity spans 10 5 in LbHa- The implied [O III] correction 
factor, using Lb i = 130 LbHa (Paper I) is noted. It changes by a factor of ~ 10 over the range of luminosities shown. The mean contribution to L[o m] from 
star formation in the host galaxy (see §5.1) is marked by 'SF'. It cannot account for the observed trend, as it contributes < 10% of L[q mi ■ 



The implied A, i/L[om] values are indicated by two hori- 
zontal lines in Fig. 4. The mean L^/L^q m ] drops from 3 000 at 
LbHa = 10 44 erg s" 1 (L bol = 10 46 1 erg s" 1 ), to 300 at L hHa = 
2 x 10 40 erg s -1 (L bol = 10 42 5 erg s -1 ). At the low luminosity 
end of the Tl sample, our mean L^/L^q 1U ] is smaller by a fac- 



tor of > 10 from the Heckman et al. result. We elaborate on this 
discrepancy in §4. 

In Figure 5, we compare L[o mj with Lyy and Lx, which are 
other measures of the bolometric luminosity. The best fit for the 
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Figure 5. The relation between L[o m] and AGN luminosity in the UV and X-ray bands, in the Tl sample. Upper limits on Lx and Luv are shown as gray dots. 
The slopes of the best fitting power laws (dashed lines) and the associated dispersions are noted. Gray solid lines depict a linear relation. (Left Panel) The UV 
band is close to the peak emission of AGN, and is therefore a good indicator of L^ \ ■ The linear relation is normalized by the L[q m] /Luv ratio in the highest 
LbHa Tl objects. The L[o m] vs - ^uv relation is non-linear, implying an increase in L(, i/£[o in] with luminosity, as shown in Fig. 4. The intrinsic relation may 
be even more non-linear, since host emission increases Luv at Luv < 10 43 erg s _1 . (Right panel) The linear relation is normalized by the L[om]/Lx ratio 
given in Heckman et al. (2005), and is consistent with the best-fit slope. 



93% of the Tl objects with detections in the UV (§2.1.1) is 

ho n,]; 42 = 0.01 L ( ™ 01 , a = 0.47 dex (5) 

As neither hom] nor Luv were used in the selection of the Tl 
sample, we derive the best fit by minimizing the 2D distances of the 
data points from the fit. The noted a is the scatter in the distance of 
log Ljo in] from the best fit relation, as in eqs. 1-3. 

The L[q in] vs. Luv slope is similar to the product of the lin- 
ear LbHa vs. Luv relation (Paper I), and the 0.72 slope of LbHa vs - 
Ljo in] (eq. 3). Note that the host contribution to the UV becomes 
non-negligible at log LbHa < 42, and can bias eq. 5. A similar best 
fit, using only Tls with log LbHa > 42, gives a slope of 0.65. Since 
the UV emission occurs close to the position of the SED peak in 
AGN (e.g. Zheng et al. 1997), eq. 5 emphasizes the non-linear re- 
lation between Ljo m] and L|, i found above. 

In eq. 4 we use the 0.6 dex scatter of log Luv around the 
Luv vs. Ljo in] relation, as an estimate for the uncertainty in deriv- 
ing L^i from L[o m] . Measurement errors are ~ 0.02 dex in Luv 
(Morrissey et al. 2007), and < 0.1 dex in hom]' as noted above. 
Therefore, the uncertainty is dominated by the physical dispersion. 
Moreover, extinction along the line of sight probably has a signifi- 
cant effect on Luv (§3.7.3 in Paper I), indicating the true dispersion 
between L[o m] and Lbd might be lower. 

A similar comparison of Ljo m] and L x , on the 43% of the Tl 
objects with ROSAT detections, gives 

ho m] ; 42 = 0.05 I$*>£ 02 , a = 0.44 dex (6) 

The coefficient and dispersion are consistent with the 
i[Oni]/^2-10 keV ratio of 0.03 and dispersion of 0.48 dex 
given in Heckman et al. (2005), found for type 1 AGN selected 
by their [O III] emission (we converted L2-iokeV to Lx using 
a x = —1.5, §2.1.2). The linear relation is not surprising, given 
the similar slopes of ho m] and L x vs. LbHa ( e 1- 3 here and eq. 
7 in Paper I). Also, the known L x « ^CqL re l at i° n ( e -8- ^ ust et 
al. 2007), is similar in slope to the Ljo m] vs - ^uv relation (eq. 
5) found here, which also points to a relation between Ljo m] 
and L x which is close to linear. Eq. 6 could be biased due to the 
low X-ray detection fraction of 43% (§2.1.2). For objects with 



FbHa > 10~ 13 ' 5 erg s _1 cm~ 2 , where the X-ray detection fraction 
is 61%, the slope is 0.91 ± 0.03. Note that previous studies of type 
1 AGN found flatter L [0 m ] V S- L x slopes: 0.7 ± 0.06 in Netzer et 
al. (2006), 0.82 ±0.04 in Panessa et al. (2006) and 0.68 ±0.2 in 
Trouille & Barger (2010). 

3.4 The correlation of Lnlr/Lblr with L/L^dd and Mbh 

In the previous sections, we find a trend of decreasing L n u a /hlla 
and L[o m] /hHa with increasing LbHa - Here we explore which of 
the parameters, LbHa> *^bh or L/L E( jd, appears to be the dominant 
parameter which drives the L NLR /L BLR trends. 

3.4.1 The Tl sample 

We bin the Tl objects based on / (= log L/L E( jd) and m (= 
log M BH /M S ) in the following manner. The objects are sorted by 
m and divided into six equal size groups. Each of these groups is 
then sorted by /, and again divided into six equal size groups. This 
ensures similar statistical errors in all bins. We repeat this process 
by sorting based on LbHa and Av. 

Figure 6 presents the derived relations of the mean values of 
AiHaAbHa and L [0 m]/LbHa as a function of L b Ha, at a fixed Av 
(left columns) or as a function of m, at a fixed / (right columns). 
Error bars denote the uncertainty in the mean. In the upper row, 
objects with strong stellar absorption are disregarded since they are 
offset to lower L n Ha/Arfta values (Fig. 3). The upper-left panel 
shows that L n Ha/^bHa decreases with LbHa, as shown in Fig. 3. 
The dependence of L nH a/ibHa on Av is weak, if at all. The up- 
per right panel shows that both m and / contribute to the trend, as 
LnHa/hna decreases with increasing m, at a fixed I, and also with 
increasing /, at a fixed m. The comparison of the upper two panels 
indicates that LbHa ( or Aral) is the main driver of the L n Ha/AjHa 
ratio, and the trends with m and / are driven by the LbHa depen- 
dence. 
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Figure 6. Mean NLR to BLR emission ratios in the Tl sample, as a function of LbHa and Av (left column), and as a function of m = logMaH/M ■) and / = log L/L^m (right column). All Tl bins in each row have 
the same number of objects, as noted in the lower left corner. Same Av or / bins are connected by solid lines, with thickness increasing with Av or decreasing / (mean Av or / noted). Error bars denote the uncertainty 
in the mean value, (a, b) Objects with strong stellar absorption (§2.4.2) are disregarded. The mean L n n a /LbHa is largely set by LbHa> and independent of Av. The LnHa/^Mta decreases both with m at a fixed / and 
with / at a fixed m, indicating that LbHa (or Lbol) is the key parameter driving the LnHa/^-bHa trend, (c, d) The mean L[o m]/£bHa is set by LbHa. and independent of Av, similar to L n Ha/ibHa- Also plotted is the PG 
sample, divided into bins with nine objects each (values from Baskin & Laor 2005a, 2005b). The PG bins are connected by dashed lines, with color and thickness matched to the Tl bins. The black line shows the 42 
Tl objects with / ~ —0.3, divided into three m bins. The PG sample shows an increase of i[om]/A>Ha with Av at a fixed LbHa, as implied by EV1 in Boroson & Green (1992), and in contrast with the Tl sample. 
On the other hand, the PG L[o mj/^bHa is consistent with Tl at m < 8.5, while the trend of decreasing L[o ni] /£bHa with increasing m reverses at m > 8.5. There are few m > 8.5 Tls, thus inducing the apparent 



discrepancy in the left panel. The reversed trend of L[o m]/£bHa is due to the dominance of radio loud AGN at m > 8.5. The la difference in L[o i 
lowest Lbiia, corresponds to the Seyfert / LINER transition at / ~ —2.7 found by Kewley et al. (2006). 
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In the lower panels, the analysis is repeated for L[o m] /^bHa- 
The behavior of L[o nil /^bHa with AGN characteristics is similar to 
the behavior of LnHa/AjHa- The trend of decreasing Lyo In ] /L(,h« 
with increasing m or / are driven by the L b Ha dependence. The 
L[q in] /i-bHa is dependent on Av only at the lowest L b H« ( see be- 
low). 

We conclude that the decrease in relative NLR luminosity in 
the Tl sample is mostly due to an increase in AGN luminosity, 
and not due to an increase in L/L^dd or A^bh- This result appar- 
ently contradicts earlier results, e.g. for the PG sample (Boroson & 
Green 1992). We therefore repeat our analysis on the PG sample, 
as described below. 



3.4.2 The PG sample 

Boroson & Green (1992) performed a PCA on the PG quasar sam- 
ple (Neugebauer et al. 1987). Their first eigenvector (EV1) indi- 
cates that L[Om]/L bH| g correlates with FWHM(H/3), while their 
EV2 indicates that £[o in] /A)H/3 anti-correlates with the absolute 
V-magnitude. In Boroson (2002), EV1 was interpreted as an indi- 
cator of L/LEdd> while EV2 follows L bo ]. As seen in Fig. 6, in the 
Tl sample we recover the EV2 trend of L[o in] /i-bHa with L bo i, hut 
do not see a trend of L[o ni] /A>Ha with Av, apparently discrepant 
withEVl. 

We add the 81 z < 0.5 PG quasars to Fig. 6, divided into 
3x3 bins and color coded as the Tl sample. We use L[om]> 
L mp , and FWHM(H/3) from Baskin & Laor (2005a, 2005b), as- 
suming FWHM(H/3)= Av and L bHa (Tl)/L bHi3 (PG) = 2.45 (cal- 
culated from common objects). The PGs have a bin with / = —0.3, 
which does not appear in Tl. Therefore, we also plot the positions 
of the 42 Tl objects with Z ~ —0.3, divided into three m bins (these 
are the 42 objects with the highest / in the I = —0.7 bin of Tl). 

In the lower left panel, it can be seen that L^q m] / ^4>Hce in- 
creases by a factor of up to ~ 2.5 over a factor of 4 in Av at a 
given L^Ha, as implied by the EV1 relations, and in contrast to the 
lack of trend in ni] l^bHa vs - Av in the Tl sample. In the lower 
right panel, it can be seen that L[o m] /A>Ha in the PG sample is 
consistent with the Tl relations at m < 8.5. However, the trend of 
decreasing L[o m] /^bHa with increasing m reverses at m > 8.5. We 
note that only 5% of the Tls have m > 8.5, compared to 33% of the 
PG quasars. 

What can cause this change of trend in L[o m]/A>Ha at high 
Mbh? A possible suspect is radio loudness, since strong radio emis- 
sion is known to be associated with enhanced narrow line emission 
(de Bruyn & Wilson 1978, and citations thereafter) and with large 
M BH (Laor 2000). Specifically, 52% of the m > 8.5 PG quasars at 
z < 0.5 (i.e. the Boroson & Green sample) are radio loud, com- 
pared to 2% of the m < 8.5 PGs (fig. 2 in Laor 2000). Therefore, 
we suspect that emission line filaments associated with radio lobes 
in radio loud quasars provide additional sites for [O III] emission at 
high m, and drives the L[o m] /AjH/J vs - Av correlation of EV1. At 
m < 8.5, the dominant trend of ni] /A)Ha m the PG sample is 
with L bo i, as found for the Tl sample. 



3.4.3 LINERs 

Kewley et al. (2006, hereafter K06), showed that SDSS type 
2s occupy two 'clouds' in the BPT diagrams (Baldwin et al. 
1981, Veilleux & Osterbrock 1987), such that LINERs (Heck- 
man 1980) have distinctly weaker L[q m] /AiH/3 an d stronger 
L[Sn]/LnHa an d i[Oi]/^nHa than Seyferts. K06 showed that 



Seyferts and LINERs are separated by a threshold of L[o mj/crf ~ 
10~ 1,9 Lq (km s~')~ 4 . This threshold is equivalent to / = —2.7, 
using the Giiltekin et al. (2009) M BH — rr* relation, and eq. 4 at 
L[0 m] = 10 6 ' 8 Lq, the typical m] at the Seyfert/LINER transi- 
tion (fig. 16 in K06). Do we see this transition in the Tl sample? 

In the lower left panel of Fig. 6, the mean Ljq m] /AjHa is 
lower by a factor of two (difference of 7rj) in the log L\^ a = 
41, logAv = 3.8 bin, compared to other bins with the same L b Ha- 
Objects in this bin have / = —2.6, the lowest / in the Tl samplfl 
and similar to the threshold I found by K06 for the Seyfert/LINER 
transition. Therefore, the drop in L[ hi] /i-bHa is probably related 
to the LINER phenomenon. The LINER interpretation is also sup- 
ported by the high L[o ij/AiHa an d £-[s n] /L n Ha ratios in the mean 
spectra of these objects (Fig. 2). Thus, the L[o m] °= L®^ relation 
found above is probably applicable only at / > —2.6, for AGN clas- 
sified as Seyferts. 



4 COMPARISON WITH PREVIOUS STUDIES 

We find here that at L bol = 10 46 erg s -1 , L bol /L [0 m] =3000; while 
at L bo i = 10 42 - 5 erg s' 1 , L hol /L [0m ] = 300 (Fig. 4). At low AGN 
luminosity, our L bo i/L[o nil differs by a factor of > 10 from results 
of previous studies, such as the L bo ]/L[o ni] — 3 500 found by Heck- 
man et al. (2004, hereafter H04). In this section, we compare our 
methods with prior work, to explore the reason for the discrepancy. 

Croom et al. (2002), H04 and Netzer et al. (2006), found that 
the mean EW([0 III]) remains constant with luminosity in type 1 
AGN samples with z < 0.3, which should overlap well the Tl sam- 
ple. At higher z and higher AGN luminosity, Sulentic et al. (2004) 
and Netzer et al. (2006) found a decrease of EW([0 III]) with L bol , 
dubbed the 'disappearing NLR'. This transition between the flat 
and decreasing parts of the EW([0 III]) vs. L 5l0O ^ relation can be 
seen in fig. 9 of Ludwig et al. (2009). 

In Paper I, we derived the mean Lhost / ^agn ratio at 5 1 00A, for 
different L bol (fig. 13 there). At log L bol = 45.5, L host / L AGN ~ 0.2, 
while at log L bo i = 43, Z. bos t / ^AGN ~ 10. Therefore, the host con- 
tribution to the continuum increases from ~ 15% to ~ 90% over a 
range of 10 2 - 5 in L bol . This increase in host contribution can also 
be seen in the decrease of stellar absorption features EW with lumi- 
nosity (fig. 8 in Croom et al 2002; fig. 18 in Paper I). Hence, if the 
host is not accounted for, at log L bo ] = 43 the measured EW([0 III]) 
will be lower by a factor of ~ 10 than its intrinsic value. Above 
log L bo ] =45.5, where the host contribution at 5100A is negligi- 
ble, the measured EW([0 III]) becomes equal to the intrinsic value. 
None of the studies mentioned above accounted for the host contri- 
bution to the continuum when measuring EW([0 III]). Therefore, at 
low luminosities (or low z), the increase in relative host contribution 
with decreasing L bo ] cancels the intrinsic increase in EW([0 III]), 
yielding the measured apparent constant EW([0 III]) vs. L bo i rela- 
tion found in earlier studies. At high L bo i or z, the host contribution 
is negligible, and the intrinsic decrease in EW([0 III]) emerged. 
Here, we use L b Ha as a bolometric indicator, which is not sub- 
ject to host contamination, and therefore the intrinsic decrease in 
EW([0 III]) is apparent at all luminosities. 



Note that due to our binning method, the lowest / bins in the upper right 
panel have a larger mean / ~ —2.2, and do not show a similar drop. 
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5 THE PHYSICAL MECHANISM BEHIND THE 

Amlr Abol VS. L bol TREND 

In §3 we have established a trend of decreasing Lnlr/^blr with 
increasing Lblr, or equivalently, a dominance of intermediate type 
AGN at low LbHa- I n this section, we explore the possible physical 
mechanisms that drive this trend. 

Changes in the ionization and density of the NLR gas will 
have a weak effect on the Ha recombination line, which is emit- 
ted proportionally to the ionization rate in the low density NLR gas 
(e.g. Osterbrock & Ferland 2006). Therefore, such changes can- 
not drive the trend in LnHa/Arol- Evidence supporting or opposing 
other possible physical mechanisms is presented in the following 
subsections. 

5.1 Star formation 

Narrow line emission can come from H II regions in the host 
galaxy, powered by ionizing radiation from young OB stars. Can a 
relative increase in star formation (SF) with decreasing L\, \ explain 
the increase in Lnlr/^blr by a factor of 10? The SF contribution 
to [O III] emission in high metallicity galaxies is relatively weak, 
and AGN hosts generally have high metallicities (Groves et al. 
2006, and references therein). Specifically, Moustakas et al. (2006) 
show that the expected L[o rn] /L n Ha from SF is 0.07, in contrast to 
the Tl sample which has L[Om]/L n Ha = 1-5 at log LbHa = 43.5, 
and L[o iu]/^ n Ha = 0.9 at log LbHa = 40. Therefore, SF does not 
dominate the narrow line emission at low AGN luminosity. The ex- 
pected contribution to L[o m j /LbHa from SF, based on the host UV 
emission in the Tl sample and the L[o ml /Luv ratio of SF galaxies, 
is shown in Fig. 4. Its derivation is given in appendix C. 

5.2 Optically thin dust 

A possible interpretation for intermediate type AGN is that the cen- 
tral source is viewed through optically thin dust close to the center, 
which extincts the AGN continuum and BLR, but has a negligible 
effect on the NLR, which originates from larger scales. Dust ex- 
tinction has been found in intermediate type AGN (Goodrich et al. 
1995, Maiolino & Reike 1995, Trippe et al. 2010), and also in the 
Tl sample (Paper I). Higher dust column densities, which obscure 
the broad optical Balmer lines, but reveal broad IR Paschen and 
Brackett lines have also been found in AGN (Veilleux et al. 1997). 
Therefore, a range of dust opacities exists in AGN. Can this range 
explain the trends seen in Figs. 3 and 4? 

We test this scenario in Figure 7, which compares Luv /LbHa 
with i-bHa/^nHa in the Tl sample. If optically thin dust drives the 
range in LbHa/L n Ha> we should see a relation between Luv /LbHa 
and -LbHa/^nHa> as most dust models have larger opacity in the 
UV than near Ha (e.g. Laor & Draine 1993). The expected slope 
of the correlation for Milky Way dust is shown in the ploQ. Clearly, 
the observed range of Ljjv/LbHa can explain only a small range in 
the observed LbHa/^nHa values. Other dust compositions, such as 
those found in the SMC and LMC, have steeper extinction curves, 
and therefore create a smaller range in LbHa/^nHa for the observed 
range in Luv /LbHa- Hence, they are even more discrepant with 

2 To calculate that effect of dust on Luv, we assume an unabsorbed UV 
slope of —0.3, apply the extinction law, redshift by the Tl sample mean z 
of 0.13, and convolve the result with the response functions of the GALEX 
filters. Our conclusions do not change for a different reasonable z or intrinsic 
UV slope. 



the data. This result coincides with the conclusion of Trippe et al. 
(2010) that most intermediate type AGN are not viewed through 
optically thin dust. 

5.3 Partial Obscuration 

Another mechanism that can reduce the observed BLR emission 
is an optically thick obscurer that partially covers the BLR. The 
following quantitative argument can be applied to decide whether 
partial obscuration is a likely mechanism behind the majority of 
intermediate types. 

A partial obscuration of the BLR in a significant number of 
objects requires an absorbers size, r a t, s , which is comparable to the 
size of the BLR, j"blr- The typical Av ~ 3000 km , implies that 
the BLR arises from a distance of /"blr ~ 10 4 GMbhA' 2 from the 
center, while the accretion disk UV continuum arises from a dis- 
tance of r uv ~ 10-100 GM BH /c 2 . Thus, r abs /ruv ~ 100- 1000, 
and the absorber will either completely absorb the UV source, or 
not absorb it at all. For example, if 90% of the BLR is covered, 
for most configurations of the obscurer we expect 90% of the ob- 
jects not to be detected in the UV. This is clearly not seen in Fig. 
7. Objects with LbHa/AiHa ~ 30 have 94% UV detections, while 
objects with LbHa/LnHa ~ 3 have 92% UV detections. Therefore, 
it is unlikely that the trend in L nH a/LbHa vs - AjHa is d ue to partial 
obscuration. 



5.4 Variability 

The BLR samples the AGN emission on timescales > 1 000 times 
shorter than the NLR (e.g. Laor 2003). Therefore, if an AGN is ob- 
served at a temporarily low state, it can appear as an intermediate- 
type AGN. However, such a scenario would imply that intermediate 
types should be relatively rare (to keep the mean luminosity high), 
contrary to the results of Figs. 3 and 4. Therefore, variability cannot 
explain the observed trend in L n Ha/^bHa vs - LbHa- 

5.5 Reddening in the NLR 

Kauffman et al. (2003) find a correlation between NLR reddening 
and AGN luminosity in type 2 AGN (fig. 21 there). They derive the 
amount of reddening from the Balmer decrement, and the AGN lu- 
minosity from L[o in] , corrected for the reddening. Could reddening 
of the NLR at high luminosity type 1 AGN, which does not affect 
the BLR, explain the drop in L„ Ha /L b Ha at high LbHa? 

If the extinction in the NLR comes from a dust screen external 
to the NLR clouds, it is unlikely that the average BLR is not af- 
fected, as the BLR originates from a region ~ 10 3 times smaller 
than the NLR, as mentioned above. Alternatively, emission line 
photons may be absorbed by dust embedded in the NLR clouds, 
which do not reside along the line of sight to type 1 AGN. How- 
ever, such dust is also unlikely to explain the trends seen in Figs. 3 
and 4, since L[o ml and L n n a , which have different dust opacities, 
show a similar trend with LbHa- 

5.6 Covering factor 

Another mechanism that can explain the decrease in the NLR to 
BLR luminosity ratio with L b Ha is a decrease in CFnlr, °r an in- 
crease in the BLR covering factor (CFblr)- In Paper I, we showed 
that the mean CFblr remains constant with L^ua, and therefore 



12 Jonathan Stern and Ari Laor 




Figure 7. The effect of reddening and NLR covering factor on Luv /i*Ha vs. LbHa/^nHa in the Tl sample. Upper limits on Luv are marked by down-pointing 
blue arrows. The symbols of other Tl objects are based on the measurement of L^wa, as described in Fig 3. The diagonal gray line depicts the effect of 
modifying the column density of dust along the line of sight to the central source (marks note steps of 0.1 in Lr-v, starting from an arbitrary value Lg_ v ). 
The horizontal line depicts the effect of modifying CFnlr. The observed range of Luv /^-Mta indicates that optically thin dust can explain only a small part of 
the range in LbHa/^-nHa- The range of LbHa/AiHa in the Tl sample is most likely dominated by a range of CFnlr • 



cannot explain the observed trend. As all other mechanisms sug- 
gested in §§5.1 - 5.5 are unlikely, we conclude that the drop in 
AiHa AbHa and L [0 m] /^bHa with L hHa is most likely due to a de- 
crease in CFnlr with AGN luminosity. A change in CFnlr will not 
affect Luv/A>Ha> as observed in the Tl sample (Fig. 7). This result 
is consistent with the result derived by Baskin & Laor (2005b) for 
the PG quasar sample, based on photoionization modeling of the 
relative strengths of [O III] A5007, [O III] A4363, and the narrow 
H/3 line. A similar suggestion was made by Ludwig et al. (2009). 
Below we derive the values of the CFnlr at different L\, \ (§6.4). 



6 IMPLICATIONS 

6.1 Intermediate type AGN 

Above we found that at L(, i < 10 44 erg s _1 , the majority of type 
1 AGN are intermediate types. The abundance of intermediate type 
AGN at low Lbol is consistent with the result of Ho et al. (1997), 
that at LbHa ~ 10 39 erg s~ 1 the LbHa is typically ~ 30% of the total 
Ha+[N II] luminosity. 

Consequently, Seyfert 1.5 samples will tend to have lower in- 
trinsic luminosities than Seyfert Is, while Seyfert 1.8s will have 
even lower L\, \. The tendency of Seyfert 1.8s to have low Z^ol has 
been noted by Deo et al. (2007, 2009) and Trippe et al. (2010). 
At luminosities below those spanned by the Tl sample, the broad 
lines may be too weak to be detectable, and unobscured AGN will 
be misclassified as obscured AGN (see next section). 

Also, we show in §5.2 that the inner region of low luminosity 
intermediate type AGN is not viewed through optically thin dust. 
Therefore, low luminosity intermediate type AGN should be gen- 
erally treated as low luminosity unobscured AGN, and not grouped 



as likely obscured type 2s. The unobscured nature of most low lu- 
minosity intermediate types has been noted previously by Barth 
(2002). 



6.2 The L[o m ] vs. LbHa relation as an indicator of true type 2 
AGN 

Various theoretical models of the BLR predict the existence of 'true 
type 2' AGN, in which the BLR does not exist (Nicastro 2000, 
Elitzur and Shlosman 2006), or is not photoionized (Laor & Davis 
2011), as opposed to standard type 2 AGN, in which the BLR is 
obscured. Candidates are often found by looking for X-ray sources 
indicating nuclear activity, with no sign of obscuration in the X-ray 
spectrum, and no sign of a BLR in their optical continuum (Rigby 
et al. 2006, Trump et al. 2009, Shi et al. 2010, Tran et al. 201 1). 

In paper I we used the L x vs. LbHa relation to test whether the 
lack of detection of broad lines is indeed significant, and justifies 
the true type 2 identification. Here we describe how the L[o m] vs. 
L[,Ha relation (eq. 3) can also be used for the same purpose. We 
compare the upper limits on L(,Ha 01 six true type 2 candidates 
with their predicted L(,Ha based on their observed L[q rn] • 

We use the three objects in Tran et al. (201 1), NGC 4594 and 
IRAS 01428-0404 from Shi et al. (2010), and Q2131-427 from 
Panessa et al. (2009), since they all have measurements of L[o mj3 
and of the continuum flux density near Ha, and estimates of Mbh 
and L/Lgjj. For each object we set an upper limit on the broad 



3 The L [0 in] of IRAS 01428-0404 is derived from the 6dF spectrum (Jones 
et al. 2004, 2009). 
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Ha flux equal to 10% of the continuum flux density near Ha 
multiplied by the expected Av. The expected Av is estimated from 



Av = 1850 (4^ 
v 10 8 M Q 



L Edd 



°- 24 km 8" 



(V) 



which is derived from eqs. 2 and 3 in Paper I. We neglect the de- 
pendence on LbH a , which has an index of 1/45. 

The six objects are plotted in Fig. 8, using the measured L[o ml 
and derived upper limit on LbHa- The Tl sample is shown in the 
background, together with eq. 3 and the associated dispersion. Ex- 
cept Q213 1-427, all candidates have upper limits which are above 
their expected LbHa- The Q213 1-427 upper limit is 2.5o~ lower than 
the expected L D Ha> but still within the distribution spanned by the 
Tl sample. Thus, the absence of a broad Ha in these six objects 
is not highly significant in any of the objects. A high angular reso- 
lution spectrum is required to exclude the strong host contribution 
near 7000A, and to be able to detect or exclude the expected weak 
Ha. Such weak Ha features have been detected from the ground 
using a 1" x 1" aperture (Barth 2002) and with HST (Ho et al. 2000, 
Barth et al. 2001). 

In Paper I, we performed a similar analysis using the Lx vs. 
L bHa relation (fig. 14 there). 1ES 1927+654 and IRAS 01428- 
0404 were found to have upper limits on LbHa which are well 
below the expected LbHa based on their Lx, suggesting they are 
potentially true type 2 AGN. These objects have high values of 
Lx/L[Ora] = 720 (1927+654) and 230 (01428-0404), which are 
3.5 and 2.4o~ above the observed mean ratio in the Tl sample 
(eq. 6). These unusually high values probably result from unusu- 
ally low L[o in] , most likely due to obscuration or absorption. This 
may therefore hint that 1ES 1927+654 and IRAS 01428-0404 are 
obscured AGN, rather than true type 2s. Obscuration of the NLR 
must be produced by galactic scale dust. Such extended dust can 
have a column of ~ 10 22 cm~ 2 , sufficient to extinct the optical line 
emission, but be completely thin to the X-ray emission near 1 keV. 

Thus, Ljo hi] ma Y provides a significant additional constraint 
on the nature of AGN, which may appear to be true type 2, based on 
their Lx and LbHa upper limit. We note in passing that the factor of 
~ 10 smaller column required to absorb the 5000 A emission, com- 
pared to the ~ 1 keV emission (e.g. Laor & Draine 1993), allows in 
principle to produce type 2 AGN where the optical BLR emission 
is absorbed, but the absorber is transparent to the X-ray emission. 
Only in the mid IR, the dust opacity becomes low enough that the 
absence of X-ray absorption implies no absorption. 



6.3 CF NLR 

In this section we derive CFnlr by assuming that in the NLR 
most incident ionizing photons are absorbed by the gas, and not 
by the embedded dust. This assumption is relaxed in the fol- 
lowing section. The low density of the NLR clouds ensures that 
an Ha recombination photon is emitted roughly 1/2.2 times per 
ionization event, independent of the NLR conditions. Hence, the 
Ha emission rate equals the emission rate of the ionizing pho- 
tons multiplied by CFnlr/2.2. The rate of ionizing photons is 
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Figure 8. The L[o in] vs. LbHa relation as a probe for true type 2 AGN. 
The Tl objects are shown as gray dots. The best power law fit and disper- 
sion (eq. 3) are indicated by black and gray lines. Also shown are upper 
limits on LbHa in six true type 2 candidates from Panessa et al. (2009), 
Shi et al. (2010), and Tran et al. (2011). These objects have an apparently 
unobscured X-ray emission, but no broad Ha detection. The upper limit 
on LbHa depends on the expected Av, which is derived from the Mbh and 
L/Lem estimates for each object (eq. 7). Except Q2131-427, all candidates 
have upper limits which are consistent with their expected LbHa based on 
IfOin]. The Q2131-427 upper limit is 2.5o~ from the expected LbHa, and 
within the distribution spanned by the Tl sample. Thus, the absence of a 
broad Ha in these six objects is not significant. 



on the mean observed vL v (1450A)/Lbea ( e( l- 5 in Paper 1), ad- 
justed to vL v (1200A) using a local slope of -0.72 (Telfer et al.), 
and multiplied by a factor of two to account for the mean amount 
of UV extinction in the Tl sample (see §3.7 in Paper I). 

The observed number of Ha photons may differ from the in- 
trinsic value due to extinction by dust. The extinction can be due 
to an external dust screen, or occur within the NLR clouds. We ne- 
glect extinction from an external dust screen, since even if all the 
UV-extincting dust affects the NLR, than L n H« is reduced merely 
by ~ 15%. Dust external to the NLR that does not affect the mean 
UV emission is unlikely, as the NLR is emitted from much larger 
scales. Within the NLR clouds, 50% of the photons directed into 
the cloud are probably subject to high optical depth and converted 
to IR radiation. The other half of emitted photons are not expected 
to be extincted, as the dust opacity at optical frequencies is sig- 
nificantly lower than at ionizing frequencies (Groves et al. 2004a). 
Therefore, we assume that the intrinsic L n H« is twice the observed 
L n Ha- 

Equating the Ha emission rate to the emission rate of the ion- 
izing photons, we get 



2Lnlia =1.3xl0 54 



Lbol CFnlr 



(8) 



et al. 2002) and vL v (1200A)/L bH a = 70. The latter ratio is based 



where we used Lboi/LbHa 
get CF NLR = 0.04 at Ltd = 10 
L b0 , = 10 42 ' 5 ergs- 1 . 



10 44 ergs-' 2.2 
130 (Paper I). From eqs. 1 and 8, we 



,45.5 



erg s and CFnlr = 0.4 at 



4 In NGC 4450 and NGC 4579, HST detected broad Ha with flux den- 
sities which are ~ 10% of the stellar continuum in 2" x 4" ground based 
measurements, where the broad Ha is invisible (Ho et al. 2000, Barth et al. 
2001). 



6.4 Why does CFnlr decrease with Lb i? 

Why does CFnlr drop with luminosity in AGN? A straightforward 
option is that the solid angle subtended by the gas surrounding the 
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nucleus decreases with increasing L\, Q \. Alternatively, the fraction 
of ionizing photons absorbed by dust embedded in the circumnu- 
clear gas may increase with L DO i, thus decreasing the measured 
CFnlr, which was derived assuming all ionizing photons are ab- 
sorbed by the gas. In this section, we quantify the ratio of dust to 
gas absorption of a given cloud, and then discuss how to differenti- 
ate between these two possible drivers of the CFnlr vs. L\, \ trend. 

6.4. 1 Dust bounded clouds 

For a given density, the gas opacity is roughly proportional to , 
where U is the ionization parameter, compared to the dust opacity 
which is independent of (J. Therefore the ionization rate, and con- 
sequently the recombination line flux, will increase with (J only up 
to a threshold Uo ~ 0.01, where the dust and gas opacities are com- 
parable (Laor & Draine 1993). At U » Uq, dust dominates the 
opacity, only ~ Uq/U of the ionizing photons are absorbed by the 
gas, and the emitted recombination line flux is only weakly depen- 
dent on U. This local effect persists when integrating over the entire 
cloud, using the U at the cloud surfacejf] (Netzer & Laor 1993), and 
also when accounting for the effect of radiation pressure on the den- 
sity profile (Dopita et al. 2002, fig. 8 there). Hence, U > Uq clouds 
are 'dust bounded'. 

Assuming a higher Lboi implies a higher mean (J, and that a 
significant number of NLR clouds are dust bounded, then at higher 
Lbol a larger fraction of ionizing photons are absorbed by dust, and 
the measured CFnlr will decrease. Evidence for the existence of 
these clouds was presented by Netzer & Laor (1993), which showed 
that dust bounded clouds naturally explain the spatial gap between 
the BLR and NLR, and why the intermediate region emits mainly 
in the IR (the 'torus'). Also, dust bounded clouds can account for 
the small scatter in U observed in the NLR (Dopita et al. 2002, 
Groves et al. 2004b). 

6.4.2 Is the CFnlr vs. Lb i trend due to dust bounded clouds or 
a drop in £2 with L DO i ? 

We denote the covering factor of the circumnuclear gas, which can 
be different from CFnlr due to dust absorption, as £2. How do we 
differentiate between an £2 independent of Lbol- coupled with dust 
bounded clouds, and a decreasing £2 with Lbol? As mentioned in 
§6.3, ~ 50% of the emission line photons are inward bound into 
the cloud, and will probably be converted to IR emission. There- 
fore, the fraction of ionizing photons initially absorbed by H I, de- 
termined by U, will have a weak effect on the final amount of IR 
energy emitted from a cloud. Specifically, the IR emission is inde- 
pendent of U up to a factor of ~ 2, in contrast with the strong effect 
of U on the emission lines described above. Hence, by comparing 
Lir/L do i with L DO l one can measure £2 directly. If dust bounded 
clouds dominate the trend seen in Fig. 3, then LrR/L DO i should be 
roughly constant with L 0Q) \, while if a decrease in £2 creates the 
trend, then we expect LrR/L DO i oc L^}' 3 . 

A trend of decreasing Lin/L DO \ with L DO i has been measured 
by Maiolino et al. (2007), Gallagher et al. (2007) and Treister et al. 
(2008). Maiolino et al. compared L(6.7jlm) with L(5100A) in 50 
type 1 AGN with 44 < logL bo i < 48.5, and found £2 oc LjjjJ 18 . Gal- 
lagher et al. compared L(l/xm— 100/xm) with L(0.1/im— 1/xm) in 
234 quasars, and found a decrease by a factor of ~ 2 in LrR/L op t 
over the luminosity range of log L DO \ = 44.5 — 47 (see their fig. 

5 Defined by neglecting the density gradient of the ionized zone. 



2). However, Gallagher et al. interpret this weak trend as due to 
increased reddening of the optical at low luminosity, and not as a 
change in £2. Treister et al. compared L(12/im) with the sum of 
optical and UV emission in 230 AGN, and found £2 oc L^ 14 at 
44 < log L bo i < 47.5. The Li R /L bo i oc L^J 015 found in these stud- 
ies is flatter than the L n Ha/Lbol x L^^ found here, and the Lb i 
measured are higher than the range of Lb D i spanned by the Tl sam- 
ple. If the trend in £2 implied by these studies indeed continues to 
the lower Lbol spanned by the Tl sample, than there seems to be a 
combined effect, both a decrease in £2 with Lb Q i, and a significant 
amount of dust bounded clouds. In a following paper, we compare 
Lir with Lb i m the 3 410 objects of the Tl sample. 

6.4.3 The type 1 / type 2 ratio 

The value of £2 decides the long standing question of the unob- 
scured AGN / obscured AGN ratio, and its dependence on AGN 
characteristics. We note that the results of a type 1 / type 2 ratio 
study based on a 2-10 keV selected AGN sample, are consistent 
with the IR studies mentioned in the previous section, at the same 
luminosity range (Hasinger 2008). Though, these results have been 
disputed by Reyes et al. (2008) and Hopkins et al. (2009). Further- 
more, the increase of Lnlr/Lblr with decreasing Lb i implies a 
stronger host dilution of the broad lines compared to the dilution 
of the narrow lines at low luminosity, which may cause misclassi- 
fication of an unobscured AGN as an obscured AGN (see §6. 1 and 
§6.2). If this effect is not accounted for, the fraction of type 2 AGN 
will seem to increase with decreasing Lb i, even if the true obscured 
fraction remains constant. 



7 CONCLUSIONS 

We analyze the dependence of the narrow Ha and [O III] A5007 
luminosities on AGN properties in 3 410 z < 0.3 type 1 AGN. We 
find the following: 

(i) The mean L [0m] and mean L nHa oc L^ a , for 10 4l) < L hHa < 
10 44.5 ergs -i 

(ii) Using the Paper I result that LbH« °= L DO i, this trend implies a de- 
crease in relative NLR luminosity with Lboi- 

(iii) The key AGN characteristic driving the relative decrease in NLR 
luminosity is L DO i, rather than Mbh or L/Lg^d- 

(iv) The simple power law dependence of the relative NLR luminosity 
on Lboi breaks at Mbh > 10 8 ' 5 M@, probably due to the dominance 
of radio loud objects, and at L/Le^ < 10 -2 ' 6 , probably due to the 
transition to LINERs. 

(v) The most likely mechanism behind this trend is a decrease in 
CFnlr with increasing Lboi- We derive CFnlr = 0.4 at Lboi = 
10 42.5 erg s -l ; and CFnlr = 0.04 at L bo i = 10 45 5 erg s _1 . 

The implications of the decrease in relative NLR luminosity 

are: 

(i) Intermediate type AGN dominate the type 1 AGN population at 
L bo i< lO^ergs- 1 . 

(ii) Intermediate type AGN at L DO i < 10 44 erg s _1 are generally not 
partially absorbed AGN. 

(iii) The implied [O III] bolometric correction factor in type 1 AGN 
changes from 3 000 at L bo i = 10 46 erg s" 1 to 300 at L bo i = 
10 42 - 5 ergs" 1 . 

(iv) The upper limit on L D Ha i n six true type 2 candidates is consis- 
tent with their expected LbHa based on their measured L[o m] . An 
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unusually high L x /i[o ra] ratio in true type 2 AGN candidates may 
indicate dust absorption of their line emission, rather than a physi- 
cal absence of the BLR. 
(v) The decrease of CFnlr with L\, Q \ may be due to a larger frac- 
tion of ionizing photons absorbed by dust within the gas, or by a 
decreasing covering factor of all circumnuclear gas. 

We thank the expert and knowledgeable referee for sugges- 
tions that significantly improved the paper, and Brent Groves for 
illuminating us on the typical EW of stellar absorption features. 
This publication makes use of data products from the SDSS project, 
funded by the Alfred P. Sloan Foundation, data from GALEX 
supported by NASA, and from the ROSAT Data Archive of the 
Max-Planck-Institut fr extraterrestrische Physik (MPE) at Garch- 
ing, Germany. 
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Appendices 



APPENDIX A: HOST SUBTRACTION 

As described in §2.2, we adjusted the Yip et al. (2004) first eigen- 
spectrum (ESI), which is used to model the host galaxy. In the 
following wavelengths, we replace the flux densities of the Yip 
et al. ESI, with the appropriate flux densities from the Hao et 
al. (2005) ESI: 4835A^t880A; 4985A-5020A; 6285A-6320A; 
6540A-6600A; 6700A-6750A. We use the 0.06 < z < 0.12 Hao 
et al. ESI, normalized by 0.014. 

In Figure Al, we show an example of the fitting procedure 
near Ha, on a host dominated spectrum. 
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Figure Al. An example of the fitting procedure near Ha, on a host domi- 
nated spectrum. The plotted spectra are the observed SDSS spectrum, the fit 
host contribution, the host-subtracted spectrum, and the emission line fits. 
Note how the prominence of the narrow Ha increases after subtracting the 
host. 



data and the fit at 5007 ± 10^0 In 99% of the objects, the total error 



on Ft 



[O in] 



is < 30% of F, 



[O m] 



Out of the remaining 43 objects, in 



32 the [O III] is not blended with other components, so we derive 
F[0 in] by simply summing the continuum subtracted flux. In four 
objects only an upper limit on F[o m] can be deduced (see below). 
In three objects only the continuum was badly fit, so we use the 
measured F[o m j . The remaining four objects are marked in Table 
1 as having a large error in F[o ni] • 

In 16 objects (0.5% of the Tl sample) the fit [O III] flux den- 
sity is < 3.5 times the local flux density error. The tabulated F[o ml 
of these objects are upper limits, derived by assuming a Gaussian 
profile with this flux density and the width of the narrow Ha. 



APPENDIX B: NARROW LINE MEASUREMENTS 

The fits of the narrow Ha and [O III] A5007 are described in §2.3 
- §2.4, and the measured fluxes are listed in Table 1. Here, we dis- 
cuss the typical errors in the measurements, the derivations of upper 
limits when the emission lines are weak, and the sensitivity of the 
fluxes to the parameters of the algorithm. 



Bl [O III] 

The [O III] is fit by a 4 th -order Gauss-Hermite (GH). In high S/N 
spectra, the profile may be too complex for the GH to achieve an 
acceptable fit. Therefore, we add to the formal error (derived from 
the spectrum errors published by SDSS) the difference between the 



B2 Narrow Ha 

The error on F n H« is calculated as described above for Fjo m0 In 
97% of the Tl sample, the relative error is < 30%. Of the 107 re- 
maining objects, in 40 only an upper limit on F n u a can be deduced. 
In three objects only the top of the broad Ha was badly fit, so we 
use the measured F n u a . The remaining 64 objects (1.9%) have poor 
fits, though due to their small number we do not attempt to further 
improve the algorithm. They are marked in Table 1. 

Our algorithm can robustly detect the narrow Ha if its flux 
density is above three times the local flux density error. The 73 

6 The mean relative offset between the data and the fit is 2%. 

7 The mean relative offset between the data and the fit is 4%. 



Type 1 low z AGN. II. NLR relative luminosity 



objects with lower flux densities are marked as upper limits, derived 
by assuming a Gaussian profile with this flux density and the width 
fit to the other narrow emission lines. 

Also, we note that the algorithm chooses the lowest x 2 from 
several fit attempts, which differ by the constraints on A3 and hi\, 
or by the number of GH coefficients used to model the broad Ha 
(§2.4). We now quantify the effect on F n n a of these different con- 
straints. The A3 and of the narrow lines are either set to zero, set 
to equal the [O III] values, or are optimized during the fit. In 68% 
of the objects the difference in the derived F n ea between the three 
options is < 0.03 dex, while in 90% of the objects the difference is 
< 0. 18 dex. Choosing a different number of GH coefficients for the 
broad Ha component changes F n u a by < 0.02 dex in 68% of the 
objects and by < 0.07 dex in 90% of the objects. 



APPENDIX C: SF CONTRIBUTION TO L [Q m ] 

In this section we evaluate the mean SF contribution to L[Om], 
in half-decade L^ a bins of the Tl sample. The mean host Ltjv 
was derived in Paper I, by subtracting the expected AGN UV 
emission, based on L\, Ha , from the mean observed L uv . To con- 
vert Luv(host) to L[q in] (host), we compare these two quantities 
in the 66 SDSS star forming galaxies described in §3.3 of Pa- 
per I. These SF galaxies have complete UV detections, and span 
10 41 - 5 < L uv < 10 44 erg s" 1 and 0.01 < z < 0.25. We find 

log L [0 m] ; 43 = 0.77 log L UV; 43 - 3 . 1 5 (CI) 

where luminosities are given in units of 10 43 erg s _1 . The implied 
mean SF contribution to L[q m] is shown in Fig. 4. 



